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Abstract 

We study the long-time behavior of variants of the telegraph process with position- 
dependent jump-rates, which result in a monotone gradient-like drift toward the ori¬ 
gin. We compute their invariant laws and obtain, via probabilistic couplings argu¬ 
ments, some quantitative estimates of the total variation distance to equilibrium. Our 
techniques extend ideas previously developed for a simplified piecewise deterministic 
Markov model of bacterial chemotaxis. 
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1 Introduction 


1.1 The model and main results 

Piecewise Deterministic Markov Process (PDMP) have been extensively studied in the 
last two decades (see [6, 7, 15] for general background) and have recently received renewed 
attention, motivated by their natural application in areas such as biology [22, 8], commu¬ 
nication networks [9] or reliability of complex systems, to name a few. Understanding the 
ergodic properties of these models, in particular the rate at which they stabilize towards 
equilibrium, has in turn increased the interest in the long-time behavior of PDMPs. 

In this paper we pursue the study of these questions on PDMP models of bacterial 
chemotaxis, initiated in [10, 11] by means of analytic tools, and deepened in [13] and [20] 
on simplified versions that can be seen as variants of Kac’s classic “telegraph process” [16]. 

We consider the simple PDMP of kinetic type (Z t ) t>0 = ((Y t ,Wt)) t>0 with values in 
M x {— 1, +1} and infinitesimal generator 

Lf(y, w) = wd y f(y, w ) + (a(y)t {yw < 0} + b(y)t {yw>0} )(f(y, - w ) - f{y , w)), (1) 

where a and b are nonnegative functions in M. That is, the continuous component Y evolves 
according to = Wt and represents the position of a particle on the real line, whereas 
the discrete component W represents the velocity of the particle and jumps between +1 
and —1, with instantaneous state-dependent rate given by a{y) (resp. b(y)) if the particle 
at position y approaches (resp. goes away from) the origin. This process describes, in a 
naive way, the motion of flagellated bacteria as a sequence of linear ’’runs”, the directions 
of which randomly change at rates that depend on the position of the bacterium. The 
emergence of macroscopical drift is expected when the response mechanism favors longer 
runs in specific directions (reflecting the propensity to move for instance towards a source 
of nutriments). We refer the reader to [22] for a scaling limit of the processes introduced 
in [10, 11] that leads to simplified models like (1). 

In the particular case where the jump rates are constants such that b > a > 0, the 
convergence to equilibrium of the process (1) has been investigated in a previous work [13], 
where fully explicit and asymptotically sharp (in the natural diffusive scaling limit of the 
process) bounds were obtained. We also refer the reader to [20] for the study of related 
models in the circle, relying on a spectral decomposition, and to [21] for a general approach 
to some kinetic models including the above one, based on functional inequalities. 

In the present work we will consider position dependent jump-rates which throughout 
will be assumed to satisfy: 

Hypothesis 1.1. Function b (resp. a) is measurable, even, non decreasing (resp. non 
increasing) on [0, +oo), bounded from below by b > 0 (resp. a> 0). Moreover we assume 
that b{y ) > a(y) for all y / 0. 

In the sequel, b stands for sup y>0 b(y) G [b, oo] and sgn : M —> {—1, +1} denotes de function 

sgn (y) = t{ y > o} - l{ y <o}- 

Let us denote by /jf’ w the law of Z t = (Y t ,Wt) when issued from Zq = (y, 
following is a our main result: 

Theorem 1.2 (Convergence to equilibrium). There exists n > 0, I\ > 0, and A 
that for any y, y G M and w, w G { — 1, +1}, 

< Ke Kly \ v \y\ e ~ xt . ( 2 ) 

TV 


y.w 

Ft 


y.w 

Mt 


w). The 
> 0 such 
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The constants above can be expressed in terms of the functions a and b following the 
lines of the proof. We will try to provide as explicit as possible bounds in each of its steps. 

The proof of Theorem 1.2 relies on a probabilistic coupling argument, reminiscent of 
Meyn-Tweedie-Foster-Lyapunov techniques, see [19, 17]. Variants of this type of methods 
have been developed in several previous works on specific instances of PDMP [4, 2, 1, 18]. 
The model under study in the present paper is harder to deal with, since the vector fields 
that drive the continuous part are not contractive. 

Our approach will be based on extensions of some ideas and methods developed in [13]. 
However, due to the non constant jump-rates we will have to work with explicit couplings 
of jump-times, for two copies of the process found at different positions. Moreover, we 
will need to make use of some discrete time Markov process embedded in their trajecto¬ 
ries (reminiscent of [3]), in order to obtain controls of the global coupling time. These 
additional technicalities prevent us from getting estimates as explicit as in [13]. 

Before delving into the proof of Theorem 1.2, we point out the explicit form of the 
equilibrium of the process ( Y , W ) and its relation to one dimensional diffusion processes 
in a convex potential. 


Proposition 1.3 (Invariant distribution). The invariant distribution of (V, W) on M X 
{—1, +1} is given by 

p(dy,dw) = ^-e~ F ^dy 0 ^(4'_i + d + ±)(dw) 

Up z 

where Cf := dy < oo and F is the convex function 

rv 

I/ 6 M 4 F{y) = / sgn {z)(b(z) — a(z)) dz. (3) 

Jo 

The domain of the Laplace transform of y is (—6 + a, b — a) x {—1, +1}. 


Example 1.4 (Laplace and Gaussian equilibria). If a and b are constant functions, then 
p{dy, dw) = l ^- 1 ZJe~ (b ~ a ^ y][ dy <g> ^(<5_i + 5 +1 )(dw). 

If a is a constant function and b is the map y H > a + |y|, then 


A i(dy, du>) = 


1 


' y2/2 dy <g> ^(<J_i + S +1 )(dw). 


Figure 1 compares in the latter case the empirical law of Y* to its invariant measure at 
increasing time instants. 


Proof of Proposition 1.3. We first note that the constant Cf is finite as soon as the func¬ 
tions a and b satisfy Hypothesis 1.1, since z e-A b(z) — a(z ) is non decreasing and positive 
on (0, oo). Furthermore, 

1 [V 

lim - / ( b(z ) — a(z)) dz = b — a. 

y->- 1- 00 y Jo 

This ensures that the Laplace transform of p, is finite on (—b + a,b — a) x {—1, +1} and 
infinite on the complement. For any function / E C 1 on M X {—1,+1} with compact 
support one has, from the definition of F, 


Lf{y , 1) + Lf(y, -1) = dy(f(y, 1) - f(y , -1)) - sgn (y)(b(y) - a(y)){f{y, 1) - f(y, -1)) 
= 9 y (f{y, 1) - f(y , -1)) - F'(x)(f(y , 1) - f{y , -1)). 
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Figure 1: Empirical distribution of Yf starting at (5, —1) for t E {2,6,10,14,18, 22} with 
a(y) = 1 and b(y) = 1 + |y|. 


An integration by parts ensures that 

J ( d y f(y , 1) - d y f(y, -l))e~ F{ - y) dy = J ( f{y , 1) - f(y , -1))F'(y)e~ F{y) dy, 

which yields f Lf(y , w) y(dy, dw) = 0. In other words, y is an invariant measure for L. Q 

The next result is independent of the previous and can be seen as a generalization 
of the renormalisation of the telegraph process studied by Kac [16] (see also [14], [13]). 
It shows that, under the suitable scaling, process (1) behaves like the diffusion processes 
expected from Proposition 1.3: 

Theorem 1.5 (Diffusive scaling). For each N > 1 let a N ,b]y ■ K. —>• M + be jump-rates 
satisfying Hypothesis 1.1 such that y i->- ajv(y) + bjy(y) is of class C 1 and 

i) a N ( 0) + b N (0) -> oo, 

ii) b N — a N —»• 2ci and —> 2 c 2 locally uniformly for some functions ci, C 2 : M —> M 

when N ^ oo. Let {Y^ N \w ( t N) )t>o denote the process driven by (1) with a — a^^b — b ]y 
and assume that Yq —> £o ' m law as N —> oo. Then, the sequence of processes 

.= ( Y ( N ) ) 

V / 7V>1 V Tjv (')/iV>l ’ 

with t n the solution of T' N (t) = |(a JV (Y ; {^ V ^) + b N (Y^^)), t n ( 0) = 0, weakly converges in 
C'([0,oo),M) when N —> oo to the solution ft of the stochastic differential equation 

dft = dB t - (sgn(£j)ci(£t) + c 2 (£t)) dt , (4) 

where (B t ) t >o is a standard Brownian motion independent from £o- 

Remark 1.6 (Diffusion in a convex potential). The drift term in (4) is odd since on and 
bN are even. Notice from point ii) above that a^(o)+b^(o) = e 2 -^ c 2 (s)ds+o(i) un if orm iy on 
compact sets as N —>• oo. Hence ojv(y) + bjv(y) —1 oo for all > 0 by i), and c^ = 0 if and 
only if a,N{y) + &jv(y) ~ on( 0) + 6^(0) for each y > 0. Thus, any diffusion with generator 
of the form \f"(y) — U'{y)f'{y) for an even convex potential U can be obtained as a limit, 
taking for instance aN{y) = «/v(0) —> oo and 6jv(y) = o/v(0) + 2 U'(y). 
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The remainder of the paper is organized as follows. In the next subsection we briefly 
recall generalities on the coupling approach to long-time convergence in total variation 
distance. We also define therein the reflected version of process (1), a detailed study 
of which is crucial for proving Theorem 1.2. Section 2 is devoted to the study of jump 
and hitting times of the latter process. A coalescent coupling for it is then constructed 
in Section 3 and the corresponding convergence estimates are established. The proof of 
Theorem 1.2 is achieved in Section 4. Finally, Theorem 1.5 is proved in Section 5. 


1.2 Preliminaries 


In the sequel we will use the notation = meaning “equal in law to” and E( A) for an 
exponential random variable with mean 1/A. 

Recall that the total variation distance between two probability measures 77 and 77 on 
a measurable space X is given by 

II 7 ? — ^IItv = i n f {lP(Ar ^ X) : X, X random variables with C(X) = 77, C(X) = 77 

If 77 and 77 are absolutely continuous with respect to a measure v with respective densities 
/ and / then 

" 1 [ f~f 


}• ( 5 ) 


I TV 


dv = 1 — / ./ A / du. 


See [17] for alternative definitions of this distance and its main properties. A pair of 
stochastic processes ( Ut , Ut)t >0 constructed in the same probability space, for which an 
almost surely finite random time T satisfying Ut+T = Ut+T for any t > 0 exists, is called 
a coalescent coupling. The random variable 


T* = 


inf jf > 0 : U t+S = Ut+s^s > oj 


is then called the coupling time. It follows in this case that 

C(Ut)-C(Ut)\\ <V{T* >t). 

A helpful notion in obtaining an effective control of the distance is stochastic domination 
(see [17] for a complete introduction). 

Definition 1.7 (Stochastic domination). Let S and T be two real random variables with 
respective cumulative distribution functions F and G. We say that S is stochastically 
smaller than T and we write S < s to. T, if F(t) > G(t ) for any i£R. 

In particular, for a couple (Ut, Ut) as above, Chernoff’s inequality yields 


C(U t )-C(Ut) < P(T > t) < E( e 


TV 


' XT ' e~ xt 


( 6 ) 


for any non-negative random variable T such that T* < s to. T, and any A > 0 in the domain 
of the Laplace transform A 1 —> K(e XT ) of T. 

We will use these ideas to obtain the exponential convergence to equilibrium for (Y, IT) 
in Theorem 1.2, and in Theorem 1.9 below for its reflected version ( X , V) which we now 
introduce. The Markov process (( Xt , Vj)) t > 0 is defined by its infinitesimal generator: 

Af(x, v ) = vd x f(x, v) + ra(x)l { ^ = _i } + b(x)t {v=1} + {z ~° } N ) (f(x, -v ) - f(x, v)), (7) 

V Mai >0}/ 

where the maps a and b satisfy Hypothesis 1.1. The term l{ x =o}(I-{:r>o}) -1 means that V 
flips from —1 to +1 as soon as X hits zero. In other words, X is reflected at zero. 
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Remark 1.8. Given a path (( Y tl Wt)) t>0 driven by (1), a path of {{X t ,Vt)) t>Q can be 
constructed taking 

X t = \Y t \, Fo = sgn(y 0 )Wo 
and defining the set of jump times of V to be 

{t > 0 : AV t ± 0} = {t > 0 : A W t + 0} U {t > 0 : Y t = 0}. 


Since W does not jump with positive probability when Y hits the origin, one can also 
construct a path of ({Yt,Wt)) t>0 from an initial value y E M and a path (( Xt,Vt)) t>0 
driven by (7): set do = 0 and (dj) i>1 for the successive hitting times of the origin and 

{Yt, Wfi) = (-iy S gn(y)(X t ,Vt) if t E [a u a i+ i). 


Let us state our results about the long time behavior of {X, V). 


Theorem 1.9. The invariant measure of {X, V) is the product measure on M + x {—1, +1} 
given by 


u{dx,dv ) 


^-e F ^ x) dx® i(5_i +S+i){dv). 


where F and Cf are given by (3). Moreover, denoting by uf ,v the law of {X t ,Vt) when 
Xq = x and Vq = v, there exists A > 0, I\ > 0 and c > 0 such that, for any x, x > 0 and 


v,v E {-!,+!}, 


x.v 

x.v 

v t 

~ v t 


TV 


< Ke~ xt e c ( xV ^\ 


( 8 ) 


Following the lines of the proof, the constants A, I\ and c can be expressed in terms of 
the jump rate functions a and b. Let us summarize some important random times involved 
in the proof of Theorem 1.9 (all related to the reflected process): 


• T( x ,v) stands for the first jump time starting at (x, v) and T(x,v) stands for its Laplace 
transform, 


• Z/ x v \ stands for the first hitting time of (0, +1) starting at (x,v), 

• T c stands for the first crossing time of the continuous components of two paths, 

• T* stands for the coupling time (i.e. from T* on the two paths are equal forever). 

Roughly, we will let evolve both paths until the first crossing time T c (which is stochasti¬ 
cally controlled by hitting times of 0) and then couple the whole processes by using explicit 
couplings of the jump-times. Notice that it is not obvious to deduce Theorem 1.2 from 
Theorem 1.9. 


2 Basic properties of the reflected process 

2.1 Distribution of the jump times 

Let us denote by Tr XiV \ the first jump time of the stochastic process {X, V) starting from 
(Ao,Lo) = (x,v) with infinitesimal generator defined by (7). This random time satisfies 

T(x, v ) = inf |t > 0 : J c(X s ) > E 
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where E is an exponential variable with unit mean and c stands for the function b when 
v = +1 and, when v = — 1 , 



( 9 ) 

( 10 ) 


The functions A 1 and B 1 are well defined since a and b are positive functions. 

Lemma 2.1 (Law of jump times). Let x 6 M+. The random variable TL, i+1 ) is absolutely 
continuous with density given by 

The random variable Tt Xj _i\ is a mixture of an absolutely continuous random variable and 
the constant variable x. Its distribution is 


e~ A{x) 5 x + a(x - t) e -W*)-A(*-t))i [0x] ( t )dt, 
where S x denotes the Dirac mass at x and dt the Lebesgue measure on M. 

Proof. We notice that Tt x+1 } is almost surely finite since / 0 °° b(x + s)ds = +oo. Let E be 
an exponential variable with unit mean and t > 0, then 

P(T (Xi+1) >t)= b(x + s)ds < E^j = e~ & b ( x + s ) ds . 

We obtain the density of Tr x + i\ by derivation. The distribution of Ti x _^ is similarly 
obtained noting that F(E > A(x)) = e~ A ^ x \ □ 

Lemma 2.2 (Laplace transform of jump times). Let x > 0 be fixed. The Laplace transform 
°fT( x ,+i) is finite on (—oo ,b). Furthermore, if X < b(x), 



< b(x) 

~ b(x) — X 
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Proof. Let x > 0 be fixed. Since b is non decreasing, we first we notice that TL + 1 ) < s to. 
E/b(x), where E is an exponential variable with mean 1. Then its Laplace transform is at 
least defined on (— 00 , b(x)) and is bounded from above by that of E/b(x ) on this interval. 

Let us now fix A < ft. Thus there exists z > 0 such that A < b(x + z). The distribution 
of TL +1 ) conditional on TL, j+1 ) > z is equal to the distribution of z + T( x+Z)+1 ) and the 
Laplace transform of T/ Xj+ i) can be split as follows 


E 




= E 




■E 




= / 6(x + t)e^ x ~ b( - x+u ^ du dt + e fo ( x - b ( x+u )'> du E e XT (x+z,+i) 


which is finite from the definition of z. 


□ 


Lemma 2.3 (Stochastic order for jump times). For 0 < x < x, we have 


-^(x,+i) —sto. T(x,+i) and T(x,—i) > s to. T(x,—i)- 
Proof. We first consider v = +1. Let E be an exponential variable with unit mean and 
T( x+1 ) := B~ 1 (E + B(x)) — x and T(x,-i) := B^i^E + B(x)) — x 


We have 



b(x + s)ds = E = 



b(x + s)ds 


with b a non-decreasing function and x > x, then clearly T( Xj+1 ) < T/ x ,+i) (inequalities 
are strict when the jump rates are strictly monotone functions). The proof is similar for 

T(x- !)• n 


We now give helpful results in order to construct a coalescent coupling of two processes 
starting from different initial data, which is in fact a generalization of a well known result 
on exponential random variables. Indeed, when a < b are positive constants, the following 
equalities in distribution hold: 


E(b) = E[a) A E(b — a), 

E(a) = E(b) + eE'{a ), 

where random variables on the left hand side are independent and e is Bernoulli with 
parameter (1 — a/b). 

Lemma 2.4 (Decomposition of jump times, part I). For x > x > 0, the identity in law 

TWd = T( Xi+1 ) A Z + 

holds, with Z + a random variable with values in (0, 00 ] independent ofT^. + i)j such that 
P (Z + > t) = exp (—B(x + t) + B(x) + B(x + t) — B{x)) for all t G [0, x). (11) 

Moreover, there exists a coupling (T/ Xj+1 ), T/ X)+1 j) such that, almost surely, 

T (x,+1) > T {x,+l) 
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and, conditionally on {T( Xj+1 n = t} 


T(x,+\) — t + £,tT(x+t,+i) 

~ c 

with T(£ +t ,+i) = Ttx+t,+i) and ft a Bernoulli r.v. independent of T( x+f)+1 ) of parameter 

b(x + t) — b(x + t ) 


Pt-= 


b(x + 1) 


e [0,1). 


We observe that if b(x) goes to +oo as x —> oo, we have Z + < +oo a.s., whereas 
Z + = +oo a.s. if b(x) = b is constant. 


Proof. Since b is a non-decreasing function, we have for any t > 0, 


B(t + x) — B(x) > B(t + x) — B(x). 

Using the representation (9) and the memoryless property of the exponential distribution, 
we thus have for all t > 0 that 


P(T (X)+1) >t) = P (E > B(x +1) - B(x)) 

= P(T( Xj+1 ) > t)F(E > B{x + t) — B(x) — B(x + t) + B(x)). 

The first statement follows. We next check that (TU+n, Ti x ,+l)) (% + i)AZ + ,%+!)), 

with (T( X)+1 ), Z + ) as before, is the required coupling. Since 

i T (x,+l) > Z +} = { T (x,+1) > T (x,+ 1 )}) ( 12 ) 

we deduce that T (X)+1) = T {Xt+1) + (T (X]+1) - T (a . >+1) ) l{T (Si+1) >T (a; , +1) }- Thus, we just need 
to check that, conditionally on {T x ,+i — 


JL 

(l{T( 2 j+ 1 )>T (a . )+ 1 ) })T(^ |+1 ) - T( X!+1 )) = ( 6 ) T(x+ t ,+l))- 

Using (12), and the expression for the density of T1 X)+1 ) together with (11), we get 


P (T(x,+i) > T( Xi+ i) + r, T( Xj+ 1 ) > s) — P (T( Xi+1 ) > Z + + r, 
b(x + t) - b(x + t) / B ( x+r+t )_ B ( x+t )) 
b(x + 1) 


Z + > s) 

e -(fl(*+t)-B(*)) 6 ^ + Q dt 


for all s, r > 0. Alternatively, 


P (T(x,+i) > T( X;+ i) + r, T( Xj+1 ) > s) — 

roo 

/ P (% +1) - T (Xi+1) > r,T (Xi+1) > r (Xi+1) |T (Xi+1) = t) e -(B(*+t)-B(*)) b ( s + 

«/ S 


Taking derivative with respect to s in the two above integrals, one concludes by comparing 
the two different obtained expressions. Q 

The function a being non-increasing, we obtain an analogous result for T( x _iy. 
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Lemma 2.5 (Decomposition of jump times, part II). For x > x > 0, the identity in law 

T(x,-l) = l) A Z- 

holds, with Z _ a random variable with values in (0,5;] independent of T( x .+i), such that 
P(Z_ > t) = exp (-A(x) + A(x — t) + A(x) — A(x — t)) for all t € [0, x). (13) 


Moreover, there exists a coupling (T^ x _^,Ti x _ i)) such that, almost surely, 

> T(x- 1 ) 


and, conditionally on {T( x ^ = t}, 

T(x,-i) =t + XtT( x -t~ i) 


with TV x _ t x) = T( x _ t +1 ) and xt a, Bernoulli r.v. independent of Tr x _ tj+ i\ of parameter 
a(x — t) — a(x — t ) 


ct t : = 


a(x — t) 


E [0,1) if t < x and a x : = 1. 


Example 2.6 (Explicit laws for jump times). In the case b(x) = b + x with b> 0 (as in 
the TCP model studied in [9, 5]), T( x , + i\ has the density 


f(x,+i)(t) = (b + x + t)e 


(t+x+b ) 2 — (x+b ) 2 


and an everywhere finite Laplace transform given by E[e AT (** + 1 l] = 1 + A rj(x + b — X), 
u 2 

with rj(u) = e” 2 " y/2n(l — $(u)) and $ the cumulative distribution function of a standard 
Gaussian variable. We also notice in this case that for 0 < x < x, 


T (x,+i) = T (x,+i) A E(x - x), 

for E(x — x) an exponential variable of mean l/(x — x) independent o/T( x +1 p and 
P(E(x -x)> T (5i+1) ) = 1 - (x - x)r]{x + b). 


2.2 Hitting time of the origin 

Let (x,v) E M+ x {—1,+1}. We notice that 

JT 

Z ix ,+i) %,_i) T B x , 

where Zr xv \ is the first hitting time of (0, +1) of a path starting from (x,v) and S x is an 
excursion above x independent of Zr x _!\. Consequently Z^ x _^ is stochastically smaller 
than Z( Xj +i). 

The Laplace transform of the hitting time of zero starting from (x,v) was explicitly 
computed in [13] in the case where a and b are both constant. Let us recall this result in 
the following proposition. 
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Proposition 2.7 (Hitting time of 0 for constant jump rates [13]). Let us define A c = 
^(Vb — \fa) 2 and, for any A < X c , 

b — a — yj (a + b — 2A ) 2 — 4 ab a + b — 2 A — yj(a + b — 2A ) 2 — 4 ab 

c( A) =-g- = -2^-‘ 

Then, for any A E (—oo, A c ], 

= e xc W and E(e xz (*’+^ = ip(\)e xcW . 

Moreover, these Laplace transforms are infinite on (A c ,oo). 

If the jump rates a and b are not constant, the evolution away from the origin is no 
longer invariant by translation. Consequently, we have to consider a new way to estimate 
the distribution of the hitting time of zero. 

Proposition 2.8 (Hitting time of 0 for general jump rates). Let M > 0 such that 

j b(M) c M(y/b(M)-^L{M)) 2 / _ £ -A(m)\ 1 

Y a(M) V J 

Then, the Laplace transform of the first hitting time Z( XiV ) of ( 0,+1) starting from (x,v) E 
M x { — 1 ,+ 1 } satisfies 


E 




(xVM)(b(M)-a(M)) 1 , - , - „ 

<Ce -a- for all A < ~(y/b(M) - y/a{M )) 2 , 


where C > 0 is an explicit constant depending on M, a and b. 

Proof. We first notice that f(x, v ) = e ax+ ^ v with a, (3 > 0, is a Lyapunov function for the 
infinitesimal generator A of ((X t ,Vt)) t>0 defined by (7). More precisely, we have 


a — 


&(*)( 


1 — e 


-2/3 


Af{x,+1) = f(x,+ 1 ) 

Af(x , — 1) = f(x, —1) —a + a(x) ^e 2/3 — 1 

for all x > 0. If we choose a, (3 > 0 and a compact set K = [0, M] x {—1,1} such that 

— a + a(M) — 1^ < 0 and a — b(M) ^1 — < 0, (14) 

by monotony of a and b there are p = p(a, fi, M , a) > 0 to be specified and c > 0 such that 

Af (x, v ) < -pf(x, v) + ct K {x, v ). 

The Laplace transform of the first hitting time of [0 ,M) x {—1,1} starting from (x,v), 
denoted t/ X jV \ := inf{t > 0 : (X t , Vf) E [0, M) x {—1,1}}, can then classically be controlled 
as follows: 

rt/\T^ x,v) 

fiXtAT^Vt^Je* = f(x,v) + / [AfiXM + pfiXs^e^ds + Nt^^ 

•J 0 

< fix, v) + N t AT(xv) 
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where (N t )t>o is a martingale with respect to the filtration generated by ((X t ,Vt))t>o- 
Taking the expectation in the previous inequality we deduce that, for x > M, 

]E[ e P r (x,„)l < e a(x-M)+fi(v+l)' (15) 

We next choose a,/3 > 0 in order to optimize p. For a,/3,M > 0 satisfying (14), we set 
p= b(M)(l — e -2 ^) — a A a — a(M)(e 2 ^ — 1) 

First we choose a > 0 such that 

2a = 6(Af)( 1 - e“ 2/3 ) + a(M)(e 2f3 - 1) 

and then 


P = \ [log ( 6 (M)) - log (o(Af))]. 

With this choice, we have 

2a = 6 (Af) - a(M) and p(M) = \ 
Condition (14) is satisfied for any M > 0 since 

—a + a(M) (e 2 ^ — lj = a — b(M) ^1 — e~ 2 ^ = — ^ 


- \/c 


(16) 


(17) 


- \/c 


We now obtain an estimate for the Laplace transform of Z^ +1 y Let M > 0 be arbitrarily 
fixed for the moment and A > 0 such that A < p(M). For x > M, we have 


Z(x,+l) — T (x,+ l) + Z{M- 1) 


(18) 


where tv X)+ 1 ) and Z^ M _y are independent. 

Let us denote by T\ and T 2 the first two jumps of V. From Figure 2.2 we see that, before 
a path starting from (M, —1) hits (0, +1), either 

1. Tf = M, in which case we have Zn ^_a = M; 

2. T\ < M and Ti > T\ , in which case we have 

Z(M-i) = 2Ti + < s to. 2M + r (Mj+ i) + Z^ M _ r); or 

3. T\ < M and T 2 < Ti, and then 

= Tj + T 2 + -Z(m—Ti+t 2j — i) <sto. 2M + Z( M _r), 

with Zr M _ 1 ) an independent hitting time of zero starting from (M, —1). 

Lemma 2.1 ensures that P(m,-i)0Tl = Af) = e~ A ^ M \ As a consequence, if is the Laplace 
transform of Z^ M _ x y one has 

tp[X) < e ~ A (M) e \M + (l- e -A(M)A e 2 AM E e \r (M ,+i) ^(A). 

For any A < p(M) we get, thanks to Holder inequality and (15), that 


E 


e \ r ( Mt+1 ) 


< E 


oP( m )t(m,+ l ) 


x A/p(M) < (M)/p(M). 
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Figure 2: The three different types of paths from (M, —1) to (0, +1). 


Thanks to (16), if M is chosen in order that 

I b(M) c M(^/b(Mj-^(Mj ) 2 A _ e -A(M)\ 1 

V a(M) V / 

then for any A < p(M), ip is finite and 

e -A{M) e \M 

^ _ e -A(M)y 2 \(M+ 0 (M)/p(M)) ■ 

Combining (18) and (15) with the previous estimate completes the proof in the main case 
x > M. If x < M then 

Z(x,-1) <sto. Z(M,~ 1) an d ^(x,+l) <sto. M + T(m,+i) + %{M- 1)j 

and one can conclude as in the previous case. □ 

3 The coupling time for the reflected process 

This section is dedicated to the construction of a coalescent coupling of two paths of the 
reflected process driven by (7) starting from two different initial conditions. 

3.1 The first crossing time 

Let us consider [x, v) and (x,v) two initial data with x < x. The first crossing time of two 
paths (X, V) and (X, V) starting respectively from (x, v ) and (x, v) is defined by 

T c = T c (x,v,x,v) = inf ji > 0 : X t = X t |. 

Since (Xt)t> o is continuous, T c is stochastically smaller than the hitting time of zero Z^ x v ^ 
of the initially upper path, whatever the joint law of the pair. The first crossing point Xt c 
is such that 

Xt c < sup X t - (x- x) < - (Z( x , v ) + x-x). (19) 

t£[0,Z M ] 1 

Notice that at time T c the two velocities are opposites. 
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3.2 A way to stick the two paths 

In what follows we assume that (Xq. Vq) = (x, +1) and (Ao,Lo) = (x,— 1) and construct 
two paths which are equal after a coalescent time T* (x). The successful coupling consists 
in producing a jump of (exactly) one of the two velocities V or V at a crossing time of their 
position components X and X. We will use the coupled jump-times studied in Lemmas 
2.4 and 2.5 in order to minimize the time required to do so, in the spirit of [13]. 

To be more precise, given x > 0 fixed, let us denote by 77+ and 77_ (respectively L- 
and L+) the first and the second inter-jump time lapses of the path starting from (x, +1) 
(resp. starting from (x, —1))- These random variables are constructed as follows. We first 
choose 77+ with distribution TL+ n and with distribution Tr x _i\ independently. We 
then define U- and L + in such a way that (77+,L+) and (77_,L_) have the laws of the 
couplings defined in Lemmas 2.4 and 2.5 respectively and that L + — 77+ and 77_ — L_ are 
independent conditionally on (77+,L_). More precisely, conditionally on 77+ and L_, we 
introduce two independent Bernoulli variables £ and x with 


m = M u+,l-) 


P(x = 1[77+, LX) 


b(x + 77+) — b{x — + 77+) 

b(x + 77+) 

a{x — L_) — a(x + 77+ — LX) 
a(x — L_) 


and 

+ 1{L_ 


=x} 


and two independent random variables L + — [7+ and U- — L_ with the same law as 
T x+ u + _l _ )+ i and respectively. Then we set 


L+ := U+ + i{L+ - U+) and U- := L_ + X (U- - L_). 


Figure 3 shows the four possible outcomes. Those where exactly one of the Bernoulli 
variables is equal to 1 allow us to stick the paths at time 17+ + L _ (i.e. on the rightmost 
corner of the rectangle): the velocities of the two paths are the same right after that 
instant, and the overshot length (beyond the rectangle’s corner) determined by the previous 
coupling is compatible with the law of the two marginal processes from that moment on 
(because of their Markov property). We then say that the coupling attempt succeeded, 
and this happens conditionally on (L7+,L_) with probability 

b(x - L_ + 77+) \ a(x - L_ + 77+) 
b(x + U + ) ) a{x — LJ) ^ 


( 20 ) 

Observe that the success or failure of the coupling attempt is determined by the 
Bernoulli random variables £ and X . If the coupling attempt fails, the two trajectories 
cross or bounce off of each other at time 77+ + L_ and by similar reason as before the 
(already determined) lengths (L+ — 77+) and (T7_ — L_) can be used to restart two (upward 
and downward) trajectories from x — L- + 77+, independently of each other conditionally 
on the past and consistently with the pathwise laws of each of the two processes. 

The coupling construction is now obvious: we repeat this scheme starting from the 
new crossing or bouncing point, until we succeed in sticking the two paths. Notice that 
this iterative algorithm is more efficient than the general procedure of the Meyn-Tweedie 


P(£ = 0, X = 1|77+, LX) + P(£ = 1, X = 0|77+, LX) = 


b(x - L_ + 77+) / a(x - L_ + 77+) 


+ 


b(x + 77+) 
6(77+) 


1 - 


a(x — LX) 


+ 1 - 


b(x + 77+) 


{L-= x}- 


-<M 
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Figure 3: Position of both paths after one step. 


method [19] since, after a fail, the two processes have already the same position (and still 
opposite velocities). 

We point out that the coupling scheme implemented for the reflected process in [13] 
in the constant jump rates case is a particular case of the above described scheme. Notice 
however that here, in general, the upper path starting from (x, +1) does not necessarily 
remain above the other path until the coupling time. We then cannot control the coupling 
time by the hitting time of 0 for the process starting at (x, +1) as it was done in [13]. On 
the other hand, contrary to the constant rates case where the coupling could only succeed 
right after the lower process hit 0, the coupling can now succeed at an arbitrary step of 
the scheme, though not with a probability bounded from below uniformly in x (this can 
be easily seen from formula (20) e.g. in the case case when a is constant and b(x ) = a + x ). 
Therefore, a new approach to estimate the coupling time is needed, which is developed in 
next subsection. 

3.3 Coupling time from a crossing point 

In this section we will use the notation P^ (resp. E x ) for the distribution (resp. the 
expectation) of a random variable associated with the coupling scheme given that the two 
copies started at position x > 0. 

We first observe that for fixed R > 0, the probability of success in one step can be 
bounded from below (considering the last term in (20) and taking expectation) uniformly 
over x E [0, R] by some number pr E (0,1) satisfying 

/•OO 

PR>e~ A{R) / b{u)e-^ h ( R+s)ds du. (21) 

Jo 

This suggests that the number of trials below a fixed height R > 0 required in order to get 
a successful coupling can be stochastically dominated by some geometric random variable. 
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Notice that we do not expect a successful coupling to occur only below level R. We will 
rather use the above remark in order to construct the scheme in such a way that the 
coalescent time will be always smaller than or equal to some real random variable that we 
can control in terms of geometric number of positive time lapses. 

First, we define a sequence of “rectangles” of potential trajectories of the two copies in 
the coupling scheme, on which the two copies will live at all times, irrespective of whether 
the coupling attempt has already been successful or not (of course once it has been so, 
their positions and velocities coincide from that moment on). More precisely, we define a 
discrete time Markov chain (<Iv) n >o starting at x by 


$0 = x > 

d> n+1 =d> n + T ( ^ +1) 


rpn-\-\ 


( 22 ) 


where conditionally on all the past up to (and including) time n, T an d ^^-i) are 
independent and respectively equal in law to T^ y + i\ and on the event {d> n = y}. 

Plainly, ($ n )n>i describes the height of the right-most corner of the n —th rectangle ob¬ 
tained by iterating the construction of Figure 3. Consider also the sequence of positive 
random variables (real time lengths) ( cr n ) n >o defined by 


a 0 = 0 , 

_ _ rp' 

0Yi+l — J- ( 


<n +1 


+ T) 


n+1 


(23) 


which give the (real) time-position of the rectangles’ right-most corners, and finally set 
£ n = (Ti , with the convention So = 0. Following Lemmas 2.4 and 2.5 and in order to 
determine at which attempt the coupling is successful, we introduce two sequences (£ n )n>i 
and (xn)n >l °f Bernoulli random variables, conditionally independent of each other given 


(<7fc, t f > / c )^,>o and such that for n > 0, on j<L n = y , 7) 


m+1 


-?n+1 _ 


_ J. rrrn , ^ —o' 

(3v,+l) _ B -1) “ S| 


^(Cn-t-l — 11 (crfc 5 ®k)k<n+l) — 


b(y + t) — b(y — s + t) 

Kv + t) 


and 


T . . a(y — s) — a(y +1 — s) „ 

^{Xn +1 = 1|(ffej ^fc)fc< n+1 ) = _ s ) + ^{s=y}- 

We also set £o = Xo = 1 for notational simplicity. Observe that (er n , <h n , £ n , Xn) n>0 is 
Markov process. We denote by (Rn) n >o the filtration it generates. We then define a 
sequence of random variables (« n ) by 


~ %„=i, X n=o} + %„=o, x „=i}- 

According to the discussion at the end of the previous subsection, the discrete time instant 
(or rectangle number) at which the coupling succeeds is 

p := inf {n > 1 : K n = 1} 

and the real time spent in order that this happens is 

p 

F* : = = ^2 &i- 

i=1 
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The trajectories of the two copies can be easily constructed from the previous objects, but 
we actually do not need to work with them. 

Let us now introduce the discrete random variable 


Pr := inf {n > 1 : K n = 1 and 4> n _i < -R}. 


Both p and pr are stopping times with respect to (En) n >o- Since p < pr a.s., we clearly 
have 

PR 

T* :=5>*. (24) 

1=1 


Our goal now is to exhibit an upper bound for the Laplace transform of the random 
time under ¥ x . We need to introduce the stopping time (with respect to (En) n >o) 


t r {x) := inf{n > 0 : G [0, R]}, 


(25) 


and the real time 


tr(x) 

^r R (x) := 

i =0 


(26) 


accumulated when the sequence ( < h n )n>o reaches [0, R] for the first time. Let ^P{ x , v ) denote 
the Laplace transform of Tt x ,v) for ( x , v ) £ X {—1, +1}. We have 


Lemma 3.1. Assume there exist positive real numbers R , A, /3 such that \ < (3, \ + f3 < b 
and <P(R'+i)(f3 + ^)^(r,_i)(A — /3) < 1. Then, tr < oo a.s. and 






gP x ^-x>R 


where r) := (^(r i+ i)(£ + A)<^ ( r_ 1 )(A - pj) 1 > 1 . 

Proof. For each x > R and /3 > A, from the stochastic monotonicity of the jump times 
(see Lemma 2.3) we get 


Ex 


e /3$i+Atri 


=e? x E 


e (/3+A)T X!+ i 


E 


e (A-/3)T x ,_! 


< e/3x T(R,+i)(P + \)(p(R-i){\ - / 5 ). 


(27) 


If (^(r i+ i)(/ 3 + A)y)(R _ 1 )(A — /3) < 1, we deduce from (27) that e /3 ‘ s>T R An+x ^ i = 1 ai p T R An is a 
positive supermartingale with respect to (3F n ) n> q, hence 


e Pxt *> R > E 


X 


g^ryjAn+A E2o A " a i 7 y T « An 


> E x [rf^ n ] . 


Letting n —oo in the last expectation we get by monotone convergence E x [i] TR ] < oo, 
hence tr < oo a.s.. Letting then n —> oo in the first expectation and using Fatou’s Lemma 
the statement follows. □ 


For each 7 > 0 and R > 0, we now set 

£r{ 7 ) : = sup E y [e 7 T ^’ + 1 )l Kl= o] • 
ye[o,i?] 
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Proposition 3.2 (Laplace transform of the coupling time starting at a crossing point). 
Assume that (R, \, (3) satisfy the conditions of Lemma 3.1 and moreover that £r(\+/3) < 1. 
Then, the Laplace transform of T * satisfies 





e l3x 1«>h ¥ , ( 0j+1) (A)^ {r _ 1 ) (A) 
1 - £ r (A + P) 


(28) 


Proof. Fix a; G M and R,X,/3 satisfying conditions of Lemma 3.1. From (24), we just 
need to estimate T/h To this aim, we consider the process (<f? n , d n ) n >o defined in terms of 
($ ra ,C7 n ) n >o in the following way: first set fo = 0 and f\ = tr(x) + 1 and for all n > 1, 


r n +1 = r n + Tr + 1, with t]j := inf{n > f n : G [0, R]} - f n . 


In other words, r n+ i is the index of the first attempt to couple the paths that follows the 
first (discrete) return time f n + of ( < Lfc)fc>o into [0,ii] after f n . Then, we set <3>o = x , 
do = 0, = d» T] , di = and 


Fn+1 

^n-t-l = ) &n-1-1 = 'y ] &i- 

i=Tn +1 


Thus, a n+ 1 is the sum of the real time 


T~n-\-l 1 



i=T n + l 


needed after a n in order to observe again a “rectangle corner” in [0, R], plus the time Cf n+1 
spent in one coupling attempt right thereafter. Then, G M+ is the position of the 

discrete chain (or rectangle corner) at the time instant <r n + 1 . Notice that for each i > 0, 
Ti is a stopping time with respect to the hltration (J r n ) n>0 and that, conditionally on 
Tr n n {<3? n = x}, (rjj, S”n) has the same law as the pair (tr(x), E T/{ ( x )) defined in (25) and 
(26). We can now write 

PR 

i=o 

where pr := inf {n > 1 : Kf n = 1} is a stopping time with respect to the filtration {Rf„) n >Q- 
We notice that pr < oo a.s. since the probability of fail in one step starting from a 
position x G [0, i?,] is uniformly bounded on [0, R] by 1 — pr. We then can write 


EJe 


ATih 


E> 

n= 1 




ri t - 

k Pr=R 


E> 

n =1 




" fT ' H 


R'T n 0,Av-f 


1-1 — — 1 


(29) 


On one hand, we have 


E x [e A ^] =E 
< E 


A^TR (-^1 &^ (Tt r( x ) + 1 


= E 


E 


e AoVjjW+lIJ , . 
c K Tfl(a:) 


Q AS Ti? (x) 


E$ 




3 Acti 


e AS rj? (a;) 


By Lemma 2.3, E y [e At71 ] = ^^(A)^,_i)(A) < ¥>(o,+i)(A)¥>( fl _i)(A) for any point y G 
[0, i?]. Thus, using also Lemma 3.1 we get 

E,[e A ^] < ^ ( o,+i)(A)^ I _i)(A)e^ 1 ->«, 
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and 


E 


e n ^-K frt =l\Gn-l 


= E 


(fin — 1 


e A ^l K , i= i 




On the other hand, we have 


E, 


Aai e /30il^ 1>i?1 ^ = o = E e AfT1+Wl i^>«l rei= o 


and for all y £ [0. /?.] 


E,, 


e A <Tl+/30ll0 1> _R][ Ki=o 


< £r(A + P)e 


0y 


Then, again from Lemma 3.1 we get 


E, 


< £ fl (A + p)e? xl *>* 


and then, for all L = 1,..., n — 1, 


E 


e ^k e ^H k >n lK _ 0 \g kl 


= Ej 


3 A a lg Wii 4> 1>R1 


<t>k -1 

< £r(A + /3)e^^ k ~ ll ^ k ~ 1>R . 
By successively conditioning in (29), we finally have 


E,[e AT "] < J2 ^(o, + i)(A)^,_i)(A)(^(A + /?))' 

n=l 

_ e ^ xtx>Rl P(o,+i) Wf(R-i) (A) 


1 — £r(\ + P) 


for parameters as required. 


□ 


Let us now verify the existence of ( R , A, (5) such that all the assumptions of Proposi¬ 
tion 3.2 hold. Notice first that for all (5 > A > 0 and each R > 0 we have 


<P(R l) (A - P) < 


+ (/3 - A)e _ (“ +/3_A ^ 
a+P -A ’ 


(30) 


thanks to the fact that > s to. E{a) A R. Since also T( K)+1 ) < s to. E{b(R)), we 

furthermore have 

V(H 1+ 1 )(A + n< 6(fl) ^ + ;a) (31) 

for all A + (3 < b(R). 

Given A > 0, we take P > A of the form P = aX for a > (or simply a > 1 if 
b = oo). Then, we have 6 > hence we find R large enough such that 

b(R)(l-e~~ aR ) > 

a — 1 
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Thanks to (30) and (31), the assumptions of Lemma 3.1 are satisfied for all A E (0, A c ), 
where 

A c := inf | A > 0 : (a + 1)A > b(R) (l - e -(s+(a-i)A)fl) _ ^±1 

Indeed, from (30) and (31), condition <pm+u(l3 + A)^(^_i)(A — j3) < 1 holds as soon as 

b{R) (a + (a - iJAe-^"- 1 ^) < (a + (a - l)A)( 6 (ii) - (a + 1)A) 

77 b{R)(a - i)A e -(“+(“-LA)R < ( a _ i)xb(R) - (a + (a - l)A)(a + 1)A 

77 (a + 1) A < b(R) (1 - e -(a+(«-i)A)/iA _ ^L 1 
\ / a — 1 

Since the previous inequality is satisfied for A = 0, by continuity we have A c > 0; the 
function of A on the r.h.s. being strictly concave, we also have A c < oo. 

Finally, notice that by Lemma 2.3 and Holder’s inequality, for any q > 1, 

£r( 7) < (1 -P*) 1 _ 1 / V[o,+i)(97)> 

with pr E (0,1) a quantity as in (21). Taking <7 = ( 7 ( 7 ) = 7 _1 , this in turn yields, for 
each fixed R > 0 , lirri sup £ 7 ^( 7 ) < (1 — pr) < 1. Therefore, there exists \' c E (0, A c ) small 

'■y—^0 

enough such that £n((a + 1)A) < 1 for all A E (0, A' c ) . 



3.4 The coupling time for the reflected process 


Let us consider two initial data (x, v ) and (x, v), with x > x. The coalescent time T*(a;, x) 
of a path (X,V) starting from (x,v) and a path starting from (x,v) is equal to 

the first crossing time T c {x,x) of both paths plus the time spent to stick them using the 
coupling described in Section 3.2. Consequently, the coupling time is stochastically smaller 
than the hitting time Zt x ^ of the origin of the upper path (X, V), plus some remainder 
term. 

For any (R, A,/3) satisfying assumptions of Proposition 3.2, the Laplace transform of 
the coupling time T*(x,x) is bounded by 


E 


3 A T,(x,x) 


< <l(o,+i)0)P(fl~i)(A) 

1 — £r(\ + f3) 


e *T c (x,x) e /3X c t Xc> R 


where the first crossing time T c {x, x) is smaller than Z(x, v ) and the first crossing point X c 
is bounded from above by ^Z^ x v \ + x — x). Consequently, 


e AT*0,z) 

/ ^(o,+i)(A)^(ij-i)(A) m _ x)/2 

e i^+P/2) Z (x,v) 

- 

- 1-S R {\ + P) 

- 


Using now Proposition 2.8 , for 0 < A < f3 satisfying conditions of Proposition 3.2 with 
A + /?/2 < |(\/6(M C ) - yja(M c )) 2 , we finally get 


E 


g AT* (x,x) 


^ rl t P(0, + l)W t P(R,-l)W 
- C 1 -£*(* + « 6 ’ ‘ 


x(b(M c ) —a(M c )) 
2 


with C is given in Proposition 2.8 and 


M c = sup 


M > 0 : 



(32) 


This concludes the proof of Theorem 1.9. 
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4 The unreflected process 


Let us construct a coalescent coupling of two unreflected processes starting from (y, w) 
and (y, w) respectively. For a given time to > 0, the coupling algorithm is the following: 

1. Define (x,v) = (|y|,iosgn(y)) and (x,v) = (|y|, w sgn(y)). 


2. Couple two reflected processes starting at (x, v) and (x, v) as in Section 3. 

3. Let them evolve until their common hitting time of 0. 

4. Construct until that time the two associated unreflected processes starting at (y, w) 


and (y,w) as explained in Remark 1.8. The algorithm stops if, when at the origin, 
the two copies have the same velocities. Otherwise, go to Step 5. 

5. Try to couple the unreflected processes starting from (0,+1) and (0,-1) before a 
fixed time to- 

6. In case of failure, return to step 1 for two initial conditions in [—to, to] x {—1,+1}. 

The only remaining task is to analyze Step 5 of this algorithm. To that end, one has 
to study the law of (Yt, Wfl) when Yo = 0 and Wo = ±1. Let us denote by (T n ) n>0 (with 
To = 0) the successive jump times of the unreflected process. The variable S n = T n — T n _i 
stands for the n th inter-jump time. In order to lighten the computation, we restrict 
ourselves to the law after 1 or 2 jumps. 

Remark 4.1 (Jump times of the unreflected process). We can explicitly compute the law of 
the jump-times of the unreflected process. For y e M, set A(y) = 4(|y|) and B(y ) = R(|y|) 
where A and B were defined on M + in Lemma 2.1. For y > 0, the law of the first jump 
time starting from (y, —1) has the density ft y ~i) given by 



a(y — t)e a<kV ift<y, 

e ~Mii)b(t — y)e~ B ( t ~ y ' ) if t > y. 


Moreover, the survival function T( y +1 )(t) := P^+i) CTl > t) is given for y > 0 by 

F {y , +1) (t) = e-^y +t ^ B ^\ 

and for y < 0 by 



For any bounded measurable function y on M x { — 1, +1}, one then has 

E(0,-1) [g(Xt, W / it)l{ri<t,T2>t}] = E (o ,-i) [g{t - 2 <Si, +l)l{Si<t, 5 i+ 5 2 >i}] 



where 
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Similarly, 


E ( 0 , + i)b(y t , w t )i T2< t,T 3 >t] = E (o,+i)b(^ - 2S2, +i)i 5 i+s 2 <t,T 3 >t] 

rt j-t—s 2 

= / / 9(t- 2s 2 ,+l)/(0 i+ i)(si)/( Sl _i)(s 2 )i ? ( Sl _ S2 , + i)(t - Si - s 2 )dsids 2 

Jo Jo 

= j g(u,+l)hi{u)du, 


where 


M«) = f Q 2 ^/(o,+i)(®i)/(«,-i)(V^)%-^,+i)(^- a i) ds i- 

Since C((Y t , Wt)[Yo = 0, Wo = w) = C((—Y t , Wt)[Yo = 0, Wo = — w), two copies as in step 
5 of the algorithm couple before time t > 0 with probability larger than 

pt 


Et = 2 / h-i(u) A /ii(it) du > 0. 


J-t 

Remark 4.2 (Explicit lower bound). ^4 lower bound of e can be derived from the fact that 
y a(y) and y 1 -^ b(y) respectively belong to [a(t),a(0)] and [6(0), 6(f)] on the interval 
[~t,t\. 


Let us now control the total duration of the algorithm. Notice that the estimates on 
the reflected process in Section 3 do no longer depend on the initial conditions, after the 
first crossing time of the reflected copies in the compact set [— R, R] x {—1,+1}. This 
implies that, after the first iteration of Step 2 in the above algorithm, the duration of each 
step can be controlled independently of the initial data, and of the previous steps. 

Moreover, the algorithm succeeds after at most a random number of iterations with 
geometric law of parameter £t 0 . Since the duration of each step in the algorithm has a 
finite exponential moment, this is thus true for the coupling time as well. The upper 
bound of Theorem 1.2 can then be deduced. As a conclusion the bounds in Theorems 1.2 
and 1.9 depends in the same way on initial data but the rate of convergence is smaller for 
the unreflected process. 

5 Diffusive scaling 

We finally prove Theorem 1.5. Omitting for a moment the sub and superscripts for nota- 
tional simplicity, and writing 

jt := Wt + re 7 (It) - 2(a(lt)l{y t wi<o} + KYt)^-{Y t w t >o}) K (Xt)W t , 

Jt := [q j s ds and Y t := Y t + n(Yt)Wt for a given positive function k of class C 1 , we see by 
Dynkin’s theorem that the processes 


M t := Y t -J t = Y t + n{Y t )W t - J t 


and 



K(Y s )ds — 2 


Y s j s ds - 2 



K'(Y s )K{Y s )W s ds 
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are local martingales with respect to the filtration generated by (Yt, Wt). 

In fact, using f(y, w) = y + n(y)w and g(y, w) = (y + n(y)w) 2 , since w 2 = 1, we have 

Lf{y , w) = w + n'(y) - 2{a(y)t {yw <o} + b(y)t {yw>0} ) K(y)w, 

Lg{y, w) = 2w(l + n'(y)w)(y + n(y)w) - J(a(y)t {yw < 0} + b(y)t {yw>0} )yK(y)w 

= 2 n(y) + 2y(w + re'(y) - 2(a(y)\ {yw < 0} + b(y)t {yw>0} )k(y)w) + 2n\y)n(y)w. 


Integrating by parts we then get that 
M 2 = Y 2 - 2 Y t J t + J 2 

rt 


= Y 2 -2 / (Y s + n(Y s )W s )j s ds -2 J s dY s + 2 J s j s ds 
Jo Jo Jo 

= N t + 2 [ K (Y s )ds -2 [ n(Y s )W s j s ds + 2 [ k' (Y a )K(Y 8 )W 8 ds 

Jo Jo Jo 


-2 [ J s dY s + 2 f Jgjsds 
Jo Jo 

rt rt 


= N t +2 f n(Y s )ds-2 f n{Y s )W s \j s - k' (Y s )\ds - 2 [ J s dM s . 
Jo Jo Jo 

Thus, noting that 


js = k'(Y s ) + Sgn(y s ) ([2 a(Y s ) K (Y s ) - 1] + 2t {YsWs>0} [1 - k(Y s )(ci(Y s ) + b(Y„))]) 


we see for k(Y s ) = ( a(Y s ) + b(Y s )) 1 that 


M t = Y t - 


Mf - 2 


rt r 


r \ *'(Ys) + m) , __, y * ( b(Y s )~a(Y s ) \ 
Jo [(a(y) + 6(y)) 2 + [ s} \a{Y s ) + b(Y s )J 

1 .... b(Y s )-a(Y s ) 


ds — 


W t 


a(Y t ) + b(Y t ) J ’ 


a(Y s ) + b(Y s ) 


+ W s sgn(y s )- 


a(Y s ) + b(Y s )) 2 


ds 


are local martingales. 

The function t n of the statement is well defined by the Cauchy-Lipschitz Theorem, 
thanks to the assumptions on the coefficients and the fact that Yt has Lipschitz trajecto¬ 
ries. Moreover, t n is strictly increasing, the coefficients a and b being positive functions. 
Recalling the dependence on N of the coefficients, and defining for each JVgN, 


A ( "> := -l 


rt %(d'Vf»ud" ) ) 


and 


0 + &iv(£s JVj ) 




+ 


sgn {6 N) )(b N {^s N) ) ~ a N (£ N) f) 


ds 


W 


(JV) 


tj v(i) 


a iv(d ^)+^Ar(d ■*) 


-0 ™<*> h JV (d" , )+6iv(d" ) ) 


.(JV), 


we see from the previous and Doob’s optional stopping theorem that the processes 






t> o 
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are local martingales with respect to the filtration of the process (£^, W^^)t> o- Defining 

for each R > 0 the stopping time = inf{f > 0 : |£ t ^| > R of |£^| > R}, we notice 
that the hypotheses imply that 


P ( sup 

N) -1 

>e 

+ P ( sup 

Wh 


/ 

V<< 


Pt N) + / sgn(£ s )ci(£ s ) + c 2 (^ s ) ds 
Jo 


> £ 


N—>oo 


-> 0 


for all e > 0 and 

E ( sup 




Jn) _ Jn) 

Si Si- 


+ E sup 


k^r 


A 1 "* - ft-’ 


AT—>00 


-> o. 


The processes ^ N \ cx[ N ^ and thus satisfy the hypotheses of Theorem 4.1 in 

[12, p. 354] (in the respective roles of the processes X n (-), A n (-) and B n (-) therein), which 
ensures that C((^ N \t > 0)) converges weakly to the unique solution of the martingale 
problem with generator given for / G (M) by 

Gf(x) := \ ~ (sgn(x)ci(x) +c 2 (x))f'(x) 

and initial law £(£o)- 
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